Brazilian Journal of Microbiology (2009) 40: 827-837 
ISSN 1517-8382 

PHYLOGENETIC ANALYSIS ON THE SOIL BACTERIA DISTRIBUTED IN KARST FOREST 

JunPei Zhou*, Ying Huang*, MingHe Mo* 

Laboratory for Conservation and Utilization of Bio-resources, Yunnan University, 650091, Kunming, P. R. China 
Submitted: July 01, 2008; Returned to authors for corrections: October 24, 2008; Approved: May 15, 2009. 

ABSTRACT 

Phylogenetic composition of bacterial community in soil of a karst forest was analyzed by culture- 
independent molecular approach. The bacterial 16S rRNA gene was amplified directly from soil DNA and 
cloned to generate a library. After screening the clone library by RFLP, 16S rRNA genes of representative 
clones were sequenced and the bacterial community was analyzed phylogenetically. The 16S rRNA gene 
inserts of 190 clones randomly selected were analyzed by RFLP and generated 126 different RFLP types. 
After sequencing, 126 non-chimeric sequences were obtained, generating 113 phylotypes. Phylogenetic 
analysis revealed that the bacteria distributed in soil of the karst forest included the members assigning 
into Proteobacteria, Acidobacteria, Planctomycetes, Chloroflexi (Green nonsulfur bacteria), 
Bacteroidetes, Verrucomicrobia, Nitrospirae, Actinobacteria (High G+C Gram-positive bacteria), 
Firmicutes (Low G+C Gram-positive bacteria) and candidate divisions (including the SPAM and GN08). 
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INTRODUCTION 

Guizhou, a province in southwest of China, is the center 
of East-Asia developing karst area with a karst area over 
5.5xl0 5 km 2 and is the largest and the most complex 
developing karst area in the world (44). Karst easily causes 
more rocks and less soil, serious soil erosion, difficult 
vegetation restoration, frequent drought, flooding disaster and 
poor moisture storage in the soil. As a result, the 
development of the regional forests is influenced and local 
farmers suffer from slow economic development (43). 

Most forests distributed on karst landform are developed 



from soluble carbonate rocks and usually characterize by 
high surface rock exposing ratio, separate soil body, shallow 
topsoil layer and rich calcium (45). Guizhou had a forest area 
about 1.6xl0 4 km 2 , with a mean forest coverage ratio of only 
11.28% (41). In some counties, the forest coverage ratio is 
less than 5%, such as Shuicheng (2.57%), Ziyun (2.46%), 
Dafang (2.40%), and Zhijin (1.14%) (41). Due to rock 
desertification, the forest and soil area in this karst province 
decrease drastically (42). Thus, how to protect the forest is 
vital for improving the local climate and facilitating the 
economical development of this region. 

For several decades, biologists have carried out their 
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investigations on the biodiversity of plant and animal in 
Guizhou karst region (45), but few attentions were focused 
on soil microorganism. Recently, bacterial diversities in soils 
of various forest ecosystems have been studied based on 
ARDRA (33), DGGE (31), PLFA (9) analyses and so on. 
Soil bacterial composition of various vegetations had been 
reported, such as oak, beech (17), pine (24), spruce (4; 12), 
broad-leaved (10) and pristine forests (6; 22). However, few 
studies were focused on the phylogenetic diversity of soil 
bacteria of karst forest. 

Soil bacteria are an essential component of the biotic 
community in natural forests and they are largely responsible 
for ecosystem function and diversity of life, but the vast 
majority of soil bacteria still remain unknown (34). 
Enrichment-based and cultural investigations on typical 
heterotrophic microbes have shown that microbes grow in 
proportion to less than 1% of total bacteria in an environment 
(2). The combination of amplification of bacterial 16S rRNA 
genes by polymerase chain reaction (PCR) and phylogenetic 
sequence analysis by restriction fragment length 
polymorphism (RFLP) has become a powerful tool to 
investigate natural bacterial communities. The aim of the 
present study was to characterize the bacterial diversity in 
soil of karst forest and present the first knowledge to 
understand bacterial composition in such environment. 

MATERIALS AND METHODS 

Site description and sample collection 

A forest predominantly formed by the plant Quercus 
glauca (25° 47' N, 107° 48' E), located in Dushan county, 
south of Guizhou province, was chosen for study. The mean 
annual temperature, precipitation and humidity in this region 
were 18.3°C, 1320 mm and 80%, respectively. On April 15, 
2006, 30 soil samples with a distance of 10 miters from each 
other were collected from 0~10 cm depth and mixed 
thoroughly as one for bacterial community analysis. The soil 
is black and calcareous, with pH 6.5. 



Soil DNA extraction, PCR amplification and cloning 

Soil DNA was extracted with a soil DNA isolation kit 
(Catalog #12800-50, MO BIO Laboratories, Inc., USA) 
following the manufacturer's instructions. Bacterial 16S 
rRNA genes were amplified by PCR using the combination 
of bacterial primer 27f (5' AGA GTT TGA TCC TGG CTC 
AG 3') and universal primer 1492r (5' GGT TAC CTT GTT 
ACG ACT T 3') (26). The PCR reaction was performed in a 
thermal cycler with the following program: preheating at 95 
°C for 2 minutes, 25 cycles at 98°C for 1 min, 50°C for 40 s, 
72°C for 2 min and a final extension of at 72°C for 10 
minutes. The amplified products were purified using an 
agarose gel DNA purification kit (No. DV805A, Takara Bio 
Inc., Otsu, Japan). Bacterial 16S rRNA gene amplicon (ca. 
1500 bases) was then excised from a 1% agarose gel and 
eluted with the same kit. Finally, the purified product was 
ligated into the pMD 18 T- vector (Takara Bio Inc., Otsu, 
Japan) and the ligation product was transformed into 
Escherichia coli DH-5a competent cells with ampicillin and 
blue/white screening following manufacturer's instructions. 

RFLP analysis 

Inserts of 16S rRNA genes from recombinant clones 
were reamplified with vector primers M13-M3 (5' 
GTAAAACGACGGCCAGT 3') and M13-RV (5' 
CAGGAAACAGCTATGAC 3'). The purified amplifications 
were subjected to restriction fragment length polymorphism 
(RFLP) by separate enzymatic digestions with Hhal (Takara 
Bio Inc., Otsu, Japan) and Mspl (Takara Bio Inc., Otsu, 
Japan) endonuc leases following the manufacturer's 
instructions, and the digested DNA fragments were 
electrophoresed in 3% agarose gels. After staining with 
ethidium bromide, the gels were photographed using an 
image-capture system UVITEC DBT-08, and scanning image 
analysis was performed manually. 

DNA sequencing and phylogenetic analysis 

One to three representative clones from each unique 
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RFLP type were selected for sequencing. The 16S rRNA 
gene inserts were sequenced using plasmid DNA as template 
and M13-20 (5' CGACGTTGTAAAACGACGGCCAGT 3') 
or M13-RV-P (5' GGAAACAGCTATGACCATGA TTAC 
3') as sequencing primer. Sequencing was done on an 
automated ABI 3730 sequencer by Beijng Genomics 
Institute. The resulting sequences (next to the primer 1492r 
and at least 600 bp) were compared with those available in 
GenBank by use of the BLAST method to determine their 
approximate phylogenetic affiliation and 16S rRNA gene 
sequence similarities (1; 14). Chimeric sequences were 
identified by use of the CHECK-CHIMERA program of the 
Ribosomal Database Project (30), and by independently 
comparing the alignments at the beginning and the end of 
each sequence and the alignments of the entire sequence. 
Sequences differing only slightly (<3%) were considered as a 
phylotype, and each phylotype was represented by a type 
sequence (19). Nucleotide sequences were initially aligned 
using Clustal X (40) and then manually adjusted. Distance 
matrices and phylogenetic trees were calculated according to 
the Kimura two-parameter model (23) and neighbor-joining 
(36) algorithms using the MEGA (version 3.1) software 
packages (25). One thousand bootstraps were performed to 
assign confidence levels to the nodes in the trees. 

The 16S rRNA gene sequences have been deposited in 
the GenBank nucleotide sequence database under Accession 
Nos. EF141940-EF142065. 

Diversity index 

Bacterial diversity was indicated using the Shannon- 
Weaver index (H') (38). 

i=n 

H'= -^pdn pi 

i=l 

TV is the number of operational taxonomic units (OTU) or 
RFLPs, and Pi is the percentage of clones of the rth OTU or 
the relative abundance of the fth RFLP. 



RESULTS 

Diversity of soil bacteria from the karst forest 

A total of 190 recombinant clones were randomly 
selected, and their 1 6S rRNA gene inserts were subjected to 
restriction endonuclease analysis (RFLP), resulting 126 
different RFLP types. After sequencing and CHECK- 
CHIMERA analysis, 126 non-chimeric sequences were 
obtained, generating 113 phylotypes. Eighty sequences 
(63.5% of total sequences) each represented a single clone 
and 35 sequences (27.8%) each represented two clones. 
Comparative analysis of the retrieved sequences showed that 
all clones belonged to the Bacteria domain. It was 
determined that most relatives of sequences (101 sequences 
representing 154 clones, 81.1% of total sequences) were 
related to sequences from environmental clones and 31 
sequences (24.6%) had relatively low levels of similarity 
(<94%) with their closest counterparts in the GenBank 
databases (Table 1, Figure 1-3). No sequence was most 
related to bacterial sequence detected in other forests or karst 
areas in the public databases except the clone FAC10 
(DQ451449) from a Taiwan forest, relative of KF092 and 
belonging to the Acidobacteria (Figure 2). 

Phylogenetic analysis on the soil bacteria from the karst 
forest 

Phylogenetic analyses placed the 113 phylotypes in the 
following 10 groups of the domain Bacteria: the 
Proteobacteria, Acidobacteria, Planctomycetes, Chloroflexi 
(Green nonsulfur bacteria), Bacteroidetes, Verrucomicrobia, 
Nitrospirae, Actinobacteria (High G+C Gram-positive 
bacteria), Firmicutes (Low G+C Gram-positive bacteria) and 
candidate divisions (including the SPAM and GN08) (Figure 
1-3). Among them, the Proteobacteria was the largest group 
including 69 clones, followed by Acidobacteria (45 clones), 
Planctomycetes (22 clones) and Actinobacteria (20 clones) 
(Table 1). 
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Table 1. Distribution of clones and RFLP types or sequences from soil of karst forest 
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1. Proteobacteria 


69 


36.4 


47 


37.2 


89-99 


1.1 Alphaproteobacteria 


19 


10.0 


11 


8.7 


91-99 


1.2 Betaproteobacteria 


17 


9.0 


12 


9.5 


92-99 


1.3 Gammaproteobacteria 


18 


9.5 


13 


10.3 


94-99 


1.4 Deltaproteobacteria 


15 


7.9 


11 


8.7 


89-99 


2. Acidobacteria 


45 


23.7 


26 


20.6 


90-98 


3. Planctomycetes 


22 


11.6 


13 


10.3 


87-96 


4. Actinobacteria 


20 


10.5 


14 


11.1 


90-99 


5. Chloroflexi 


5 


2.6 


4 


3.2 


85-93 


6. Bacteroidetes 


7 


3.7 


6 


4.8 


94-98 


7. Verrucomicrobia 


6 


3.2 


6 


4.8 


91-97 


8. Nitrospirae 


7 


3.5 


5 


4.0 


95-98 


9. Firmicutes 


2 


1.1 


2 


1.6 


95-97 


10. Candidate division 


7 


3.7 


3 


2.4 


89-99 



Closest relatives as determined by the BLAST method. Each sequence representing a single RFLP type. S: The sequence 



similarity to its closest relatives. 



Proteobacteria 

Sixty-nine clones, represented by 47 sequences and 
accounting for 36.4% of the clone library, were 
phylogenetically associated with the following taxa of 
Proteobacteria with similarities between 89%-99%: the 
Alphaproteobacteria (number of sequences, ns=ll, number 
of clones, nc=19), Betaproteobacteria (ns=12, nc=17), 
Gammaproteobacteria (ns=13, nc=18) and 

Deltaproteobacteria (ns=ll, nc=15) (Table 1, Figure 1). 

A total of sixteen clones, represented by 10 sequences, 
were related to cultured members and belonged to putatively 
Sphingomonadales (sequence KF003), Rhizobiales (KF091 
and KF002), Rhodospirillales (KF082), Nitrosomonadales 
(KF088), Burkholderiales (KF052 and KF081), 
Xanthomonadales (KF024) and Pseudomonadales (KF046 
and KF062). Sequences KF040, KF125 and their closest 
relatives (Entotheonella clones) were grouped in a clade with 
a bootstrap confidence value of 99% (Figure 1). 
Entotheonella, a new Deltaproteobacteria genus, had been 
found thus far only in sponges of the family Theonellidae 
(37). 



Acidobacteria 

Forty-five clones, represented by 26 sequences and 
accounting for 23.7% of the clone library, were clustered 
with the uncultivated bacterial sequences of the 
Acidobacteria with similarities between 90%-98% (Table 1, 
Figure 2). Figure 2 showed a phylogenetic tree of 
Acidobacteria, which was grouped into at least 5 
acidobacterial clusters. Acidobacteria form a newly devised 
division of Bacteria, probably as diverse as Proteobacteria or 
gram-positive bacteria. The definition of this phylum was 
based on the analysis of 16S rRNA gene sequences retrieved 
from cloned rRNA genes and phylogenetically related to 
several cultivated species such as the Fe (Ill)-reducing 
Geothrix fermentans (29; 35). The clone F08 WMSP2 
(DQ450696), relative of KF099, was detected in Alpine 
tundra wet meadow soil. The clone DA038 (AJ000986), 
relative of KF049 respectively, was found from a grassland 
soil in the Netherlands. 

Planctomycetes 

Thirteen sequences, representing 22 clones (11.6% of 
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94 r KF039 |EF142017| 

Betaproteobacteria clone: 027 |AB252914| 
L Betaproteobacteria clone Ellin6095 |AY234747| 
n KF074 |EF142052| 
KF065 |EF142043| 

Hydrogenopbilaceae clone ID33 |AY578170| 
KF052 |EF142030| 
99 L Massitia timonae sp. |U54470| 
KF083 |EF142061| 
KF042 |EF142020| 
Betaproteobacteria clone HAuD-MB/2-1 |AB176695| 
Betaproteobacteria clone AKYG1851 |AY921995| 

r KF088 |EF141940| 
99 1 Nitrosomonas sp. NmS6 |AY123798| 
KF124 (KF014) |EF141976| 
KF025 |EF142003| 

Uncultured bacterium clone CR99-7-72 |AF429106| 
- Betaproteobacteria clone AKYH490 |AY921 821 1 
r KF081 |EF142059| 
"99 1 Leptothrix sp S1 1 |DQ241397| 

KF084 |EF142062| 

99 p KF095 | EF141947| 

58J L Gammaproteobacteria clone YU21-B |AB174846| 

61 P |— KF050 (KF067) |EF142028| 

H 95^ Gammaproteobacteria clone LiUU-3-334 |AY509440 

! r- KF097 |EF141949| 

9T- Gammaproteobacteria clone AKYG629 |AY921 861 1 
_ 83^KF117(KF103, KF108, KF122) |EF141969| 
50 ' " 



Gammaproteobacteria clone MERTZ2CM300 |AF424122| 

KF111 |EF141963| 

97 , — KF047 |EF142025| 
19 1 — Gammaproteobacteria clone Dover96 |AY499680 

1 Alcanivorax sp. 521-1 |DQ659430| 

^7 99 p KF046 |EF142024| 

I Pseudomonas strain 520-20 |AJ292426| 

53] p KF062 |EF142040| 

99 L CelMbrio sp. R-20952 |AJ786793| 

I Xanthomonadaceae clone SG-3 |AF548381 

1 KF024 |EF142002| 

KF070 |EF142048| - 




99 , KF079 |EF142057| 

' Alpbaproteobacteria clone IricESI |AJ566640| 

- KF055 |EF142033| 

Afphaproteobacteria clone CRIB-04 |DQ123621| 
RNzobiales isolate CSBT10/2 |AJ810382| 

KF091 (KF002) |EF141943| 
Pedomicrobium sp. |X97691 

KF113 |EF141965| 
99 L Alpbaproteobacteria clone KY40 |AB116424| 
99, KF082 |EF142060| 

' — Amaricocctis macauensis sp. |U88042 
KF075 |EF142053| 

Alpbaproteobacteria clone BICI23 |AJ3181 1 1 1 

i KF003 |EF141981| 

98 Kaistobacter terrae sp. |AB258386| 

99iKF071 |EF142049| 

^Alpbaproteobacteria Ellin361 clone |AF498743| 
KF015 |EF141993| 



9^ Alpbaproteobacteria clone AKYH1 192 |AY921758| ' 

98 r Ebtotheonella clone Dd-Ent-79 |AY897121| , 

Qq ®~>\ 1 — Entotbeonella clone Dd-Ent-69 |AY897120| 

1 KF040 |EF142018| 

KF125 |EF141977| 

KF010 |EF141988| 

— Deitaproteobacteria clone QLBB01 1 |AY862012| 
99 , KF041 |EF142019| 

' Deitaproteobacteria clone JG34-KF-383 |AJ532716 
rKF053|EF142031| 

~9§- Deitaproteobacteria clone AKYH1030 |AY921695| 

- KF005 |EF141983| 

- Deitaproteobacteria clone MND4 |AF293008| 
- KF105 |EF141957| 

- Deitaproteobacteria clone HPDOMI1F10 |AY851907| 
75, — Deitaproteobacteria clone HAVOmatl 20 |EF032748| 

-KF009 |EF141987| 
Deitaproteobacteria clone |690311| 

KF106 |EF141958| 
KF093 |EF141945| 

Deitaproteobacteria clone: Y219 |AB1 16507| 




KF033 |EF142011| 



- Tbermosipbo sp. MV1063 |AJ419874| 



69, — KF007 (KF048) |EF141 985| 

!- Acidobacteria clone AKYG951 |AY922080| 
70^ KF032 |EF142010| 

Acidobacteria clone AKYH536 |AY922135| 

H 1 KF085 |EF142063| 

' Acidobacteriaceae clone DA023 |Y07586| 

92 r KF030 (KF056) |EF142008| 
f^- Acidobacteria clone 354D |AY571790| 
5?i — KF035 |EF142013| 
67^- Acidobacteria clone AKYH1068 |AY922094| 
i KF060 |EF142038| 

100 Acidobacteria clone BAC-14A1 |AY214904| 
10Q- KF031 |EF142009| 

I L Acidobacteria clone AKYH730 |AY922152| 

i| r— KF008 |EF1419S6| 

99^ Acidobacteria clone 354C |AY571797| 

KF080 |EF142058| 

9^KF126 |EF141978| 

Acidobacteriaceae clone VHS-B3-48 |DQ394942| 
KF057 |EF142035| 

Acidobacteriaceae clone AT-s2 |AY225640| 
KF120 |EF141972| 

KF023 |EF142001| 
KF099 |EF141951| 

Acidobacteria clone F08_WMSP2 |DQ450696| 
-KF123 |EF141975| 

„i i KF004 |EF141982| 

jSn — Acidobacteria clone AKYG742 |AY922059| 
61^ Acidobacteria clone AKYG500 |AY922037| 

100 - Acidobacteria clone AKYG1109 |AY921847| 
L KF110 |EF141962| 
KF092 |EF141944| 

Acidobacteria clone FAC10 |DQ451449| 
KF049 |EF142027| 
90' Acidobacteria clone DA038 |AJ000986| 

- KF101 |EF141953| 

— Acidobacteria clone AKYG433 |AY92201 1 1 
Acidobacteria clone AKYG1139 |AY921868| 

KF061 |EF142039| 

Acidobacteriaceae clone AT-S3-24 |AY225641| 

KF109 |EF141961| 




KF089 |EF141941| 

- Acidobacteria clone: Y72 |AB1 1 6490| 

Uncultured bacterium clone CV76 |DQ499315| 

100 KF107 |EF141959| 

jOtP KF034 |EF14201 2| 



10O 



- Planctomycetes clone BBBB-1 10 |AY226234| 
100— KF018 |EF141996| 

' — Planctomyces sp |X81954| 

97 , KF051 |EF142029| 

" I ' Planctomycetales clone |AF424502| 

' Planctomycetales clone Dover382 | AY4998 1 1 1 

_^ KF096 |EF141948| 

Planctomycetales clone HPDOMI2E12 |AY851885| 
— KF090 |EF141942| 
KF038 |EF142016| 

Planctomycetales clone: K31 |AB116406| 
— KF086 |EF142064| 

Pirellula staleyi strain ATCC 35122 |AF399914| 

100 , KF016 |EF141994| 

^ Planctomycetales clone Ellin7244 |AY673410| 

KF026 |EF142004| 

KF077 |EF142055| 

KF102 |EF141954| 
Planctomycetales YNPRH54A |AF465657| 

Planctomycetales clone AKYH674 |AY921762| 

KF058 (KF069) |EF142036| 

99^ — Planctomycetales clone CJuql4 |AF239693| 

Thermosipho sp. MV1063 |AJ419874| 



87[ 




Figure 1. 16S rRNA gene-based dendrogram showing 
phylogenetic relationships of bacterial phylotypes from karst 
forest soil (shown in bold) to members of the Proteobacteria 
from public database. Bootstrap values (n=1000 replicates) of 
>50% are reported as percentages. The scale bar represents 
the number of changes per nucleotide position. Thermosipho 
sp. MV1063 (AJ419874) was used as outgroup. Sequences of 
RFLP types differing only slightly (<3%) are shown in 
parentheses. Accession numbers are given at the end of each 
sequence. 



Figure 2. 16S rRNA gene-based dendrogram showing 
phylogenetic relationships of bacterial phylotypes from karst 
forest soil (shown in bold) to members of the Acidobacteria 
and Planctomycetes from public database. Bootstrap values 
(n=1000 replicates) of >50% are reported as percentages. 
The scale bar represents the number of changes per 
nucleotide position. Thermosipho sp. MV1063 (AJ419874) 
was used as outgroup. Sequences of RFLP types differing 
only slightly (<3%) are shown in parentheses. Accession 
numbers are given at the end of each sequence. 
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Paenibacillus sp. 324 |AFD21924| 

— KF020 |EF141998| 
- Bacillus sp YY|AF414443| 

- KF054 |EF142032| 

VK, Actinobacteria clone AKYG965 |AY922081| 

97I L KF022 |EF142000| 

^ I KF104|EF141956| 

id i Actinobacteria clone LV57-1-22 |AY642555| 

\- Actinobacteria clone FBP234 |AY250866| 
1 — KF021 |EF141999| 

"-Actinobacteria clone WCHB1-08 |AF050573| 
KF044 |EF142022| 
- Kineosporia strain 1-132 |AB003935| 
KF043|EF142021| 

- Actlnobacteria clone FBP218 |AY250865| 

- Frankia sp. strain FMr16 |AY502033| 
KF119 |EF141971| 

-KF121 |EF141973| 

KF098 |EF141950| 

100 r Rubrobacteridae clone Ellin6048 |AY234700| 
1 KF066 |EF142044| 

— Rubrobacteridae clone #0649-1G9 |AF2341 19| 

— KF017 (KF012) |EF141995| 
KF068 |EF142046| 

i Actlnobacteria clone #0319-7L14 |AF234156| 

rj KF028 |EF142006| 

— KF076 |EF142054| 




100 KF094 |EF141946| 

' — soil bacterium clone S134 |AF507726| 
r KF116|EF141968| 
^jd pKF112(KF115, KF087) |EF141964| 

jgr Nitrospira strain GC86 |Y14644| 

97 Nitrospira clone JG36-GS-2 |AJ582041 1 

10[ 83^ KF027 |EF142005) 



KF073 |EF142051| 

— Bacteroidetes clone AKYG1686 |AY921971| 
r KF037 |EF142015| 
10S- Bacferojdefes clone AKYG1020 |AY921733| 
__r- KF059 |EF142037| 

100 — Pedobacter strain DSM 14825T |AJ438170| 




in AS 1.3801 DQ021 903 



KF036|EF142014| 

Flavobacterlum saliodium strs 
KF045 |EF142023| 

10$ Flavobacterlum sp |M62797| 

KF013|EF141991| 

Verrucomicrobia clone C06_WMSP1 |DQ450785| 
Spartobscteria clone Gsoil 144 |AB245342| 

- KF006 [KF019) |EF1419S4| 
KF029 |EF142007| 

KF001 |EF141979| 

Verrucomicrobla clone YNPRH34A |AF465651| 
KF072 |EF142050| 

Verrucomicrobiales clone Ellin518 |AY960781| 

Chloroflexi bacterium H5 |AF234688| 

KF114|EF141966| 
- Chloroflexi clone Nubeena191 |AY500105| 

Chloroflexi bacterium SJA-35 |AJ009460| 
10rJ KF064 |EF142042| 

— KF063 |EF142041| 
Kouleothrix aurantiaca strain MYSI-A |AB079639| 

10(J KF078 |EF142056| 



candidate division GN08 DQ329588|- 
KF118 |EF141970| 



T-KF011 |EF141989| 

100 — candidate division SPAM clone AKYG1047 |AY921949| 

r KF100 |EF141952| 

100- Uncultured bacterium clone: MIZ17 |AB179508| 

Thermosipho sp MV1063 |AJ419874| 



Figure 3. 16S rRNA gene-based dendrogram showing 
phylogenetic relationships of bacterial phylotypes from karst 
forest soil (shown in bold) to members of the Firmicutes, 
Actlnobacteria, Nitrospirae, Bacteroidetes, 

Verrucomicrobia, Chloroflexi and Candidate divisions from 
public database. Bootstrap values (n=1000 replicates) of 
>50% are reported as percentages. The scale bar represents 
the number of changes per nucleotide position. Thermosipho 
sp. MV1063 (AJ4 19874) was used as outgroup. Sequences of 
RFLP types differing only slightly (<3%) are shown in 
parentheses. Accession numbers are given at the end of each 
sequence. 



the clone library), were grouped into the Planctomycetes. 
Seven sequences were related with relatively low similarities 
between 87%-93% to cultured or uncultured bacterial 
sequences listed in the GeneBank database (Table 1, Figure 
2). Molecular microbial ecology has provided new evidence 
showing that Planctomycetes bacteria are ubiquitous and 
constitute a representative part of the natural bacterial 
population (20). The sequence HPDOMI2E12 (AY851885), 
which was clustered with KF096, was a member from 
stromatolites of Hamelin Pool in Shark Bay, Western 
Australia. 

Actinobacteria 

Fourteen sequences, representing 20 clones (10.5% of the 
clone library), were clustered with the Actinobacteria (Table 
1, Figure 3). These bacterial clones were related with 
similarities between 90%-99% to cultured or uncultured 
bacterial sequences listed in the GeneBank database. Five 
sequences were related to classified members and belonged 
to putatively the Actinomycetales (KF043) and 
Rubrobacterales (KF012, KF017, KF066 and KF068). 

Chloroflexi (Green nonsulfur bacteria) 

Four sequences representing 5 clones were related to 
Chloroflexi sequences listed in the GeneBank database (85%- 
93% similarity) (Table 1, Figure 3). The low similarities to 
these members indicated that the corresponding bacteria 
detected in the soil of karst forest belonged to putatively new 
taxonomic groups. 

Bacteroidetes, Verrucomicrobia, Nitrospirae and 
Firmicutes 

Six sequences were clustered with the Bacteroidetes 
(Table 1, Figure 3). The closest relatives of KF059 and 
KF036 (KF045) were strains belonging to the 
Flavobacteriaceae and Sphingobacteriales, respectively. 
Some species in the Flavobacteriaceae degrade soluble 
cellulose derivatives, and nine cellulolytic isolates were 
assigned to Flavobacterium johnsoniae (27). Six sequences 
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were grouped with uncultured bacterial sequences of the 
Verrucomicrobia. Furthermore, five and two sequences were 
closely related to the Nitrospirae and Firmicutes, 
respectively. KF020 and KF054 were related to classified 
members and belonged to putatively the Bacillales. 

DISCUSSION 

To date, few molecular microbiological data on bacterial 
community in soil of karst forest were presented. We firstly 
presented the knowledge on bacterial community and 
revealed a considerable number of novel and unknown 
bacterial sequences and a high diversity of putative bacterial 
communities in this environment. Sequences with low 
similarities (<94%) to bacterial sequences listed in the 
GeneBank database were mainly distributed in 
Deltaproteobacteria, Acidobacteria, Planctomycetes, 
Actinobacteria and Chloroflexi (Table 1). Clones of the 
Proteobacteria, Acidobacteria, Planctomycetes and 
Actinobacteria accounted for more than 80% of the clone 
library, while that of other 6 groups each accounted for less 
than 4%. That confirmed Proteobacteria, Acidobacteria, 
Planctomycetes and Actinobacteria were absolutely the 
dominating communities in soil of karst forest. Besides, 
Proteobacteria and Acidobacteria were two largest groups 
(more than 60% of the clone library) (Table 1). 

Diversity of bacterial communities in soil of karst forest 
was quite different with that in other forests reported until 
now (Table 2). The Shannon-Weaver index (FT) of karst 
forest (2.32) was higher than others forests, including the 
British Columbia forest. This reflected bacterial diversity in 
the ecosystem is higher than that in other forests. Karst forest 
was a relatively fragile ecosystem type (45). The high 
bacterial diversity in soil of karst forest could possibly 
contribute to its tolerance for natural disturbance and self- 
recovery and regeneration processes. 

Composition of bacterial communities in soil of karst 
forest was also unique (Table 2), which was putatively 
related to unknown soil environmental variables. Although 



the forests with broad-leaved in Ailaoshan, Xishuangbanna 
(10) and this study site are located in the Southwest of China, 
their compositions of bacterial communities are totally quite 
different (except the group of Betaproteobacteria, Table 2). 
Similar results were obtained from the spruce forests in 
British Columbia, Canada (4, 12) (especially the difference of 
Acidobacteria, Actinobacteria and Verrucomicrobia, Table 2) 
and from the pristine forests in Brazil (6, 22) (especially the 
difference of Proteobacteria, Acidobacteria, Firmicutes and 
Fibrobacter, Table 2). That seemed compositions of bacterial 
communities from soil of forest were largely independent of 
vegetation types and geographic distance. But Fierer and 
Jackson (15) discussed that bacterial diversity and 
community composition could largely be explained by soil 
pH, but unrelated to plant diversity and geographic distance, 
respectively. 

The closest relatives of KF088 and KF112, with a high 
similarity value of more than 97%, were strains belonging to 
Nitrosomonas and Nitrospira respectively (Figure 1, 3). 
Nitrosomonas is a genus of ammonia-oxidizing 
proteobacteria (39). Nitrospira strains are the dominant nitrite 
oxidizers in most environmental samples tested so far (8). 
The Rhizobiales is famous for nitrogen-fixing bacteria that 
form nodules on host plants (7). Some strains of the 
Burkholderiales are capable of nitrogen fixation, such as 
Ralstonia taiwanensis (1 1). And the closest relative of KF062 
was Cellvibrio bacteria that are capable of cellulose- 
degrading, such as Cellvibrio fulvus (5) (Figure 1). Based on 
culture-dependent method, Long et al. (28) reported 
ammonifiers (7.6><10 6 bacteria per gram soil), nitrobacteria 
(1.14*10 3 ), nitrogen fixing bacteria (1.28><10 3 ) and cellulose 
decomposers (5.68xl0 4 ) in soil of Maolan karst forest. 

The closest relatives of KF091 and KF081 were 
Pedomicrobium and Leptothrix strains respectively (Figure 
1). Pedomicrobium plays an important role in iron- and 
manganese-oxidization (16). Leptothrix is known to be 
capable of oxidizing both iron (II) and manganese (II) to 
ferric hydroxide and manganese oxide (13). Relatives of 
KF113, KF033, KF089 and KF038 are clones detected in 
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coastal marine sediment beneath areas of intensive shellfish 
aquaculture where sulfur cycle was accelerated (3). It has 
been reported that sulfur-oxidizing bacteria play a role in the 
dissolution of limestone (32). 

The Acidobacteria and the Planctomycetes were poorly 
studied phylogenetic divisions that have been detected in 
many clonal analyses and are thought to be of great 
ecological significance to many ecosystems (18; 29). But as 
limited cultivated species, their ecological functions and 
possible impacts on soils remain unclear at present. Two 
novel genera of the Planctomycetes, Candidatus Kuenenia 
stuttgartiensis and Candidatus Brocadia annamoxidans are 
capable of catalyzing the anaerobic oxidation of ammonium 
(21). 

Our present study has given a useful insight into bacterial 
populations in soil of karst forest and can be used as a 
starting point for in-depth studies. Although studies based on 
culture-independent methods make more difficult valid 
statements about the ecological role that clones might play in 
the environment, such studies will help understand the 
diversity and structure of bacterial communities in soil of 
karst forest and contributions of bacteria to the ecosystem. 
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